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Abstract. This tutorial introduces the dynamics of charged particles in a 
radiofrequcncy trap in a very general manner to point out the differences between 
the dynamics in a quadrupole and in a multipole trap. When dense samples are 
trapped, the dynamics is modified by the Coulomb repulsion between ions. To take 
into account this repulsion, we propose to use a method, originally developed for 
particles in Penning trap, that model the ion cloud as a cold fluid. This method can 
not reproduce the organisation of cold clouds as crystals but it allows one to scale the 
size of large samples with the trapping parameters and the number of trapped ions, 
for different linear geometries of trap. 

PACS numbers: 52.27.Jt, 52.25.Kn, 37.10.Pq 
1. Introduction 

Quadrupole radiofrequency traps appear to be the obvious choice for most applications 
where charged particles need to be trapped, cooled and manipulated by laser 
interactions. Optical frequency metrology and quantum information processing give 
examples of such a control of the ion motion by atom-laser interaction as shown in several 
contributions to this special issue. Up to now, all the protocols make use of the harmonic 
nature of the potential induced by the radiofrequency field and were demonstrated with 
one or few ions. On the other hand, higher order radiofrequency traps are efficiently 
used with large samples where buffer gas cooling is preferred to laser cooling. This is 
the case, for example, for microwave frequency standards working with large samples of 
Hg+ ions |1],|2] or with experiments studying cold reactions with molecules [3|,|lj. In the 
prospect of new experiments where one could think of applying quantum information 
processing methods to a larger sample (of atomic or molecular ions), it may be useful 
to consider using higher order radiofrequency traps than the usual quadrupole trap. 
The first purpose of this tutorial is to give a general background about the dynamics 
of ions in radiofrequency traps to point out the differences between the dynamics inside 
a quadrupole and a multipole trap (here we call multipole trap a linear radiofrequency 
trap with a higher symmetry order than a quadrupole). The second objective of this 
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tutorial is to present a mean-field approach already used for particles in Penning traps 
and called the nonneutral plasma. Using this method, it is possible to scale the size of 
a trapped sample with the relevant physical parameters of the problem. This general 
approach is useful to understand what is so special about the quadrupolar geometry even 
when large samples are trapped. This tutorial is organised in two parts, the first one 
describes the dynamics of single charged particle in a radiofrequency field of very general 
geometry. The second part introduces the mean field approach that allows one to take 
into account the Coulomb repulsion in large samples. This model is used in the linear 
geometry to establish relations between the trapping parameters and the characteristics 
of the trapped sample. 

2. General description of radiofrequency trapping 

2.1. The most general case 

Let's try to get some insights about the dynamics of trapped ions (mass m and charge 
q) starting with a non specific geometry. The development presented here is greatly 
inspired from a broad review written by D. Gerlich about multipoles and their use to 
study ion-molecule reaction dynamics p]. Let's assume that the potential between a set 
of electrodes can be written as the sum of a static voltage $s(r) and a radiofrequency 
voltage $r/(r, t) oscillating at frequency VL and responsible for the effective trapping of 
the charged particle. The motion of a single ion is then governed by the local electric 
field 

E(r, t) = V^sij) + V<l>./(r, t) = Es(r) + Eo(r) cos((]t). (1) 

To compute the trajectory of this ion, the difficulty lies in the integration of the equations 
of motion in a time varying electric field. So we first assume that there is no static electric 
field (Es(r) = 0). If the electric field were homogeneous, the motion would revert to the 
rf-driven oscillation 

rm = cos(fit) = -acos(()t), (2) 

where a is the amplitude of this oscillation. As Eo(r) is non homogeneous, this amplitude 
varies in space and the motion can be split into two contributions: the motion driven by 
the rf field Ri(i:) = — a(t) cos(f2t) and a slower motion Ro(t) induced by the variation of 
the amplitude of the radiofrequency field. The local electric field can then be expanded 

as 

Eo(r(t)) = Eo(Ro) - (a(t).V)Eo(Ro) cos(l]t) + ■ ■ • (3) 
The resolution proposed here for the dynamical equations relies on two approximations 

(i) the first order expansion of Eq around Rq in Eq. ([3]). 

(ii) the adiabatic approximation which assumes that the typical evolution time scale of 
a and Rq is far longer than the rf period. 
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In this frame of approximations, the slow dynamics is determined by 

^Ro = -^grad(Eo^) (4) 

and, as far as Ro(t) is concerned, it is as if the particle is trapped in a static potential 
well, called the ]5sew(io-potential: 

V'(,) = pl^, (5) 

Building a pseudopotential well implies grad(£'Q(r)) > in the three directions. In 
practice, the fulfilment of this condition is not sufficient to assure the stability of the ion's 
trajectory. A conservative stability criterion is given by a local adiabaticity criterion, 
introduced by Teloy and Gerlich [Sj [5] , and defined by the relative variation of Eq seen 
over the amplitude of the driven motion oscillation: 

|2(a.V)Eo| 2g|VEo| 
^"'^ = lEol = 

This adiabaticity parameter depends on the location in the trap except for the particular 
case of quadratic radiofrequency voltage for which VEq is homogeneous. In the 
multipole case, numerical simulations and experimental observations mentioned by D. 
Gerlich [5] show that the empirical limit rjad < 0.3 guarantees adiabaticity in most 
cases. A more recent experimental study of the loss mechanism in a 22-pole trap [H] 
complementary to a more general model of effective trapping volume for multipole [7] 
has demonstrated stability up to rjad < 0.36 ± 0.02. A safe criterion for stability could 
then be r]ad{Rmax) < 0.34. 

If the static electric field is switched on again, according to the superposition 
theorem, the motion is governed by the effective static potential 

The rf-driven motion (also called micromotion) Ri(t) is deduced from the motion in the 
pseudopotential (called macromotion): 

RiW = -^^^cos(r]t) (8) 

and its amplitude is proportionnal to the local rf-electric field. As this electric field must 
increase with Rq for effective trapping, the micromotion amplitude increases with the 
distance from the trap center. As this driven motion can not be cooled, the only way 
to reduce this amplitude is to keep the ion at the node of the electric field. 



2.2. Focus on the linear trap 

In most experiments requiring a precise control of the dynamics and kinetic energy of the 
ions, micromotion is a side effect that can bring a large contribution to the Doppler effect 
and reduces the precision reached on the kinetic energy of the ions. In a linear geometry 
where a radiofrequency voltage is applied only in the transverse plane with a translation 
symmetry along the remaining axis, the node of the radiofrequency electric field is the 
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symmetry axis. When large samples are trapped, it gives rise to a reduced micromotion 
compared to spherical traps, as the ions are closer to the node of the electric field. In 
this section, we assume such a linear geometry where the trap consist of 2k equally 
spaced rods as electrodes for transverse trapping and where the confinement along the 
axis of symmetry is reached by a static potential (the geometry of the electrodes used 
for the axial confinement is not very relevant in this tutorial). 

As a first step, we ignore this static potential and assume that voltages +V(t)/2 = 
Us/2 — Vo/2cos(f2t) and —V{t)/2 are applied on alternate electrodes (so the voltage 
difference between neighbouring electrodes is V{t), notice that in Gerlich's paper this 
voltage is chosen to be 2V{t)). The 2k electrodes are rods located at the same distance ro 
from the trap axis and solutions of the Laplace equation A$ = are linear combinations 
of 

<l>k{r,t) = %it){r/ro)'cos{ke) (9) 

where (r, 6) are the polar coordinates in the (x, y) plane. The contribution of each 
order k is fixed by the boundary conditions defined by the shape and position of the 
electrodes. Figure [1] shows the equipotential lines in the transverse plane of an ideal 
multipole where only one contribution of defined symmetry order defines the potential. 
This potential distribution can be realised ideally if the electrode surfaces exactly match 
the equipotential surfaces. Deriving the radiofrequency electric field generated by such 
a potential, the equations governing the motion in the plane of a 2fc-pole are 

x/ro - Fkit){r/rof-^cos{k -1)9 = (10) 
y/ro + Fkit)ir/ro)'-'sm{k-l)e =0 

with Fk{t) = kqUs/{2mrQ) — A;gVo/(2mrQ) cos(fit). As an example, let's have a look at 
the equations of motion of an ion in an octupole {2k = 8): 

x + F4{t){x^ -3y^x)/rl = 

y-F,{t){y^-3x'y)/rl=0. (11) 
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These equations are non-linear and coupled. The non-linearity makes the stability of 
individual trajectories sensitive to their initial conditions. As a consequence, it is not 
possible to define absolute stability conditions based only on the working parameters of 
the trap. An estimation of the local stability is provided by the adiabaticity criterion 
(Eq. ([6])). For a 2fc-pole, it increases like r'^"^ and is uniform only in the quadrupolar 
case {2k = 4): 

Vad = Hk-l)^^^. (12) 

In a quadrupole, and only in a quadrupole, the equations of motion take the simple 
forms: 

X - F2{t)x = 

y + F2{t)y=0 (13) 

which are linear and uncoupled for x and y. The introduction of the reduced time scale 
^ = Qt/2 leads to the well-known Mathieu equation 

d^ii 

— (a„ - 2g„ cos(fit))M(0 = (14) 

{u = x,y), which belongs to the family of differential equations with periodic coefficients. 
Solutions of this equation can be found in many textbooks [HI |9] . The solutions have 
the form 

u{t) = Ae*^"* C2ne^"™ + Be-''^-' ^ Ca^e"^"™ (15) 

n n 

with Uu = Pu depending only on the Mathieu parameters and through 

continuous fractions, and qu depend on the trapping parameters by: 
2gVo 2qVo 



4gf/. 4qUs 



[17) 



These solutions are stable if < /5„ < 1. This condition defines areas in the plane 
defined by (a„,g^i): the stability regions. In practice and because of technical limitations, 
radiofrequency traps are operated in the lowest stability region (the one including (0,0)). 

The solutions of the Mathieu equation can be expanded in the lowest order 
approximation which requires (|a«|,g'^) ^ 1 and implies that Pu ^ 1- In this hmit 



Pu = \lau + ql/2 and 

u{t) = U cos{uJut) (^-'- + ^ cos^lt 

This is exactly what we would obtain using the pseudopotential well resulting from 
the adiabatic approximation. It says that the main motion (or macromotion) is a 
harmonic oscillation of frequency Ux, perturbed by the driven radiofrequency motion 
(or micromotion). The amplitude of this micromotion is proportional to the harmonic 
oscillation amplitude and to the Mathieu equation parameter g^. 
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The Mathieu equations, or more generally the equations of motion given by Eq 
ffTOj) . are strictly true for a single ion in an ideal trap and do not include the Coulomb 
repulsion between ions. When several ions are trapped, the correlations between ions 
(often called space charge effect in its mean field approach) can be ignored for warm 
enough or very dilute samples where the averaged Coulomb repulsion is low compared 
to the kinetic energy. The coupling parameter T = q"^ /{ineoakBT), the ratio of the 
average nearest-neighbor Coulomb repulsion energy and the thermal energy, quantifies 
this competition between correlation and thermal motion (the usual definition sets a 
as the Wigner-Seitz radius and is related to the density by 4-7^0^/3 = 1). For low F 
(F ^ 1), the Coulomb repulsion can be ignored and it is relevant to use the single ion 
equations of motion to learn about the ions dynamics. With laser cooling techniques, F 
as high as few hundred can be reached. As soon as F > 1, correlations can not be ignored 
to calculate dynamics or equilibrium properties of the trapped ions. This problem is 
widely addressed by two methods. The first one is based on computer simulations of 
the experiment by calculating the position and velocity of each ion by a Monte-Carlo 
method [10] , when only the statistical equilibrium is studied, or by molecular-dynamics 
simulations pTl [T2| [T3] when information about the dynamics is also wanted. These 
computer simulations are an efficient tool to calculate the crystal structure appearing in 
trapped ions as soon as F > 100, depending on the trapping parameter, the number of 
ions and the temperature [121 [lH [iS, [IE] . If is even possible to include buffer gas or laser 
Doppler cooling in the model, by an effective damping force [13|, [T7] or by a stochastic 
force induced by the momentum kicks in each collision or absorption/emission process 
^18j . In this tutorial I wish to present another method which has been developed for 
electrons or ions in Penning traps but which is very relevant for ions in radiofrequency 
trap, using a model for the atomic sample called the nonneutral plasma. This model uses 
a mean-field approach where the ion-ion Coulomb repulsion is taken into account by the 
total field created by the charge distribution. One can show [I9] that the equilibrium 
state of such a system in a harmonic pseudopotential is equivalent to that of a One 
Component Plasma (OCP), a model system where the charged particles are embedded 
in a neutralizing background charge. It allows one to calculate the density and size of 
trapped charged samples in the dense and cold limit, a requirement which fits very well 
with laser cooled atoms. Contrary to the molecular dynamics or Monte Carlo methods, 
it can not reproduce the internal Coulomb crystal structure observed in experiments 
but it is of great help to scale the size of a linear trap with the desired number of 
ions. This is the question that I try to answer in the following section. Most of the 
concepts introduced here can be found in the review by D.H.E Dubin and T.M. O'Neil 
[T9] concerning the thermal equilibrium states of trapped nonneutral plasmas. 

3. The mean- field approach 

When a large sample is trapped, a mean field approach can be used to study the global 
behaviour of the sample. This treatment has been mostly used in the context of Penning 



About the dynamics and thermodynamics of trapped ions 



7 



traps, for electrons or ions [201 EI] but also in 3D radiofrequency quadrupole traps (or 
Paul traps) [22] . In this model, a thermal equilibrium state is assumed and the sample is 
treated like a cold fluid. This approximation holds if the plasma size is large compared 
to the Debye length X^,. This length characterises the scale over which a mean field 
approach is relevant to treat the Coulomb repulsion. It depends on the density n and 



temperature T of the plasma as = yksTeo/q'^n and decreases for cold and dense 
plasmas. The cold fluid model allows one also to calculate the aspect ratio (radius over 
length) of a cloud trapped in a harmonic potential [23] . These predictions have been 
confirmed experimentally with laser cooled ions in a Penning trap [23] and in a linear 
rf quadrupole in the isotropic [25] and anisotropic [26] case. 

In the context of radiofrequency trap, the equilibrium state concept requires us 
to use the pseudopotential approach to represent the trapping potential. Like in the 
previous section, we assume a linear geometry but the same points developed here could 
be applied to 3D multipole trap. To find an experimental characterisation of the thermal 
equilibrium of an atomic cloud trapped in a 3D octopole trap, the reader should refer to 
[27] where the density profile is studied by means of laser induced fluorescence resolved 
in space. In the linear geometry, the pseudopotential in the transverse plane of an ideal 
multipole is defined by (see Eq (jTj)) 

rW = ^f n"" + ?^n'cos*». (19) 



32Sk Vro/ 2 Vro. 

where Sk = mVL^rl/ {2k'^) is a characteristic energy. In the following, we consider that 
there is no static voltage Ug applied to the rods but we take into account the static 
potential required for axial confinement and its effect on the transverse pseudopotential. 
The total effective potential becomes 



y*i^-) = \^[-] +^(2.^-0, (20) 



where k is a loss factor depending on the geometry of the end electrodes relative to 
the rods. It includes all screening effects that can explain the reduction between the 
potential applied on the end electrodes Vend and the one effectively seen by the ions. For 
2k = 4, the shape of the transverse pseudopotential remains quadratic and no major 
impact is expected from the transverse deconfining effect of the axial confinement. On 
the contrary, for 2fc = 8, 12 . . ., the shape of the pseudopotential is modified. The axis of 
the trap becomes an unstable position and the potential minimum is shifted to r = r^m 
defined by 

2k-4: _ '^0 ^^^kl'^Vend 

Wik-l)qVo' 

We see later that this shift of the potential minimum from the center of the trap has no 
effect on the density distribution calculated in the mean-field model. 

According to the work of Dubin, DriscoU, O'Neil and Prasad [2Sl [221 US] and 
assuming the ergodic hypothesis, the thermal equilibrium state of a nonneutral plasma 
with non negligible correlations can be described by the iV-particle Gibbs distribution 



2fc-4 _ '0 ^'J'-'ki^^ end /^-rN 

min 2 ! h. 1'\„T/2' \ ) 
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n(r) = A/'exp 



(22) 



for plasmas as small as > 100 (for smaller N, fluctuations are too big to identify the 
average over the microcanonical and the canonical-or Gibbs-distribution) . Integrated 
over the N velocities and — 1 positions, the Gibbs distribution gives access to the 
plasma density distribution n(r): 

kBT_ 

where A/" is a normalisation constant, and £{y) is the energy of a particle at position r. 
It results from the contribution of the confining potential ^*(r) = g07-(r) and the mean- 
field Coulomb repulsion potential q4>q{jc) created by all the charges surrounding point r 
(the image charge effect is neglected as we assume the charges are far enough from any 
conducting surface [30]). The normalisation issue can be solved by defining the density 
with respect to the density in the center of the trap n(0) = no : 

nir) = noexp[\E'(r)] (23) 

^(r) = - ^ (0T(r) + 0,(r) - 0^(0) - 0,(0)) . (24) 

The density distribution and the potential created by the charges are related by the 
Poisson equation 

AMr) = (25) 

This system of equations is then self-consistent but one can get rid of some normalisation 
difficulties by studying only the logarithmic density profile \E'(r) [28l [29] and using 
Eq. ( I2^1l25i) . one can show that 



exp[^(r)] - —A(Pt 
quo 



(26) 



It is obvious from Eq. fl26l) that the Debye length relative to the central density 



= ykBTeo/{q'^no) will be the relevant length scale for the density profile. This 
equation also shows that the static contribution to the effective potential V* has no 
impact on the density profile as this contribution obeys the Laplace equation A0 = 
and that only the pseudopotential associated with the radiofrequency field controls this 
profile. Eq. (!26l) also points out the difference between quadrupoles and multipoles, 
regarding the density profile of a nonneutral plasma. Indeed, in the cold fluid limit 

(T 0), to prevent divergence of the density, exp[\E'(r)] ^AcpT must also tend 

to which means that n(r) —>■ eoA0T/g. This equation can also be deduced from 
the condition for mechanical equilibrium of a shell of ions for which the trapping field 
balances the field created by the charges inside the shell. In our radiofrequency multipole 
context, this becomes 

lim n(r) = ' ° — (27) 



which is uniform only for the quadrupolar geometry. For higher order geometry, the 
density is expected to increase with the distance from the center of the trap, leading to 
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an empty center and a geometry like a tube. This is confirmed by molecular dynamics 
simulations presented in [31] for calcium ions in an octopole. 

We now want to use the differential equation ( 126|) to scale a cloud size. It appears 
that the situation is very different for a quadrupole from other multipole geometries. 
So we first deal with quadrupole geometry, which is very close to the Penning trap 
configuration, and then expand our method to higher order geometry. 



3.1. In a linear quadrupole 

In the particular case of a linear quadrupole, the pseudopotential can be written as 



^^(r, z) = -uy + -uy (28) 

with u"^ = — Ux being defined by the solution of the Mathieu equation [8]. As 

mentioned above, the contribution determining the density profile is 

AMr,z)=2-u;l (29) 

and does not depend on the static end voltage. The second term in the differential 
equation (!26|) is then a constant 2meQUjl./ {q^uo) that can be understood as the ratio of 
the uniform limit density Uc = 2meoUjl/q'^ that is reached for low temperature, divided 
by the realised central density uq. This ratio is bigger than 1 as the limit density Uc is 
the highest reachable density for a given harmonic potential. In the following this ratio 
Uc/uq is expressed as 7 + 1, 7 measuring how far the system is from the uniform limit 
density (7 > 0). Finally, the differential equation determining the density distribution 
is 

A^(r, z) = ^ [exp[^(r, z)] - 7 - 1] (30) 

To simplify the integration, we introduce the reduced coordinates (p, = / Xd, z/ Xd) 
and solve: 

AM/(p,0=exp[M/(p,0]-7-l- (31) 

Case of a prolate cloud L ^ R: In the following, we assume the cloud is prolate 
and the dependence of the density with z is negligible compared to its dependence with 
r. This simplifies the problem to the integration of a one dimension differential equation 

1 9^ d^^ 

To have an insight into the role of 7 in this equation, it is useful to look for an 
approximate solution, valid close to the trap center. For \I'(p) ^ 1, exp[\l/(p)] — 1 ~ \E'(p). 
In this approximation, the solution of the differential equation (1521) is 7(1 — /o(p)) where 
/o(p) is the modified Bessel function of order 0. Using the lowest order expansion of 
/o(p), one can deduce that for p^ -C 4/7, n{p) ~ no{l — 7p^/4). So the smaller 7 is, the 
larger (in reduced parameter) is the range where the density profile is nearly fiat. 

Integration of Eq (132|) (with the limit conditions \l/(0) = and \I/'(0) = 0) results 
in the profile ^I/(p) and its exponential n{p)/no which implicitly depends on 7. To get 
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real values for the number of trapped ions, the absolute density profile n(r) and the 
size of the sample, one has to find other relations between the relevant functions and 
parameters: 

• First, in the frame of our approximation, the total number of ions is 



J J n(r)2Tcrdrdz = 2L J n(r)2TTrdr (33) 



y"exp[^(p)]27rpdp. (34) 



if L is the half length of the cloud. By using the reduced parameter p and the 
definition of A^i, one can show that 

iV ^ keTeo 
2L g2 

This number of ions depends implicitly on the central density through the 
parameter 7 = Uc/uq — 1. Here, the number of ions per unit length N/2L, the 
temperature T and the central density uq are the three relevant parameters for the 
density profile problem and only two of them are free as they are related by eq . 

Second, to relate the two physical free parameters N/2L and T to the real size of 
the cloud, it is useful to notice that, even if Ad does not appear explicitly in the 
differential equation, it is not a free parameter since 



7 + 1 2mc(j2 



(35) 



X 



As an experimentalist, one may want to know what the size of an ion cloud will 
be. The solution is there: for a given temperature T and a given number of ions per 
unit length {N/2L), 7 has to be found to fulfil the condition fixed by equation (IMll . 
Then is deduced from equation fl35l) and the radial size of the cloud is given by 
R = XoPmax- The strength of the harmonic pseudopotential \fmh)x plays a role only in 
this scaling. It becomes necessary to know its value only when the size of the cloud is 
required. Examples of integrated radial profiles in the reduced and full scale coordinate 
are shown on figure [2] for the same number of ions per unit length and same potential 
well vTiiJl. but for different temperatures. This figure illustrates the differences between 
warm and dilute cloud where pmax — 1 (T = 10000 K, case a and b) and a colder and 
denser cloud where Pmax ^ 1 (T = 5 K, case e and f), where the cold fluid model is 
relevant. 

As mentioned above, theoretically, it is not appropriate to model warm and dilute 
samples by the cold fluid model, nevertheless, we can use the radial size calculated in 
this condition as an indication and compare it to the radius of the same sample cooled 
down. This is very important for experimental issue, as the actual number of laser cooled 
trapped ions can be limited by the size it takes just after or during the ionisation process. 
Indeed, when ions are created by electron bombardement, their initial temperature can 
reach an order of magnitude of 10 000 K and comparison of plots (b,d,f) of figure [2] shows 
that the radial size of the same sample at 10 OOOK is nearly 4 times the size at 300 K 
and 6 times the one at 5 K. This illustrates one of the advantage of photoionisation over 
electron bombardement to ionise a neutral beam. During the photoionisation [32], the 
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20 40 12 3 



P r(m) X10-' 

Figure 2. Density profiles versus the reduced radius p (a,c,e) or the radius r (b,d,f) 
for a prolate cloud with 10^ ions per mm, at different temperatures (a,b: T=10 000 
K, 7 = 5, c,d: T=300 K, 7 = 0.04, e,f: T=5 K, 7 = 10-^^). The scaling factor Xd 
between p and r is calculated for Ca"'' ions in a trap with uJx/^tt = 1 MHz. Notice that 
the X-axis of figure b has a different scale from d and f, and that the x-axis of figure 
e has a different scale from a and c . 



ions can be continuously loaded and laser cooled and the sample never reach as high 
temperature as by electron bombardement. Then the cloud remains smaller and more 
ions should be loaded in a trap, with less radiofrequency heating [33l [M] . 

For a laser cooled nonneutral plasma {N/2L > 10^/mm, T < 1 K), or for a room 
temperature dense sample {N/2L > 10^/mm, T = 300 K), 7 is negligible in Eq (1351) 
which can be simplified to 



X 



This shows that in the cold fiuid limit, for a given set (T, N/2L), the size of a cylindrical 
cloud scales as (ma;^)~^/^. 

For temperatures close to the Doppler limit, the density profile of large enough 
sample (> 10^ ions/mm) is fiat and expands to values of p ^ 1. We can then 
approximate the density profile by exp[\E'(p)] = 1 for p < p^ax and exp[^'(p)] = 
for p > pmax- The relation between the number of ions and the temperature is then 
simplified to: 

-^Pmax- (37) 



2L q 



2 
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As a consequence, the radial size R of the cloud is given by: 

1 ./IE (38) 

which is independent of the temperature and can be considered as the minimum limit 
radius Rm for a sample. This limit for the size of a prolate cold cloud is equivalent to a 
maximum limit for the density: 

which is of course the limit density Uc we introduced at the beginning of the text. 

To have an idea of this size, for a laser cooled nonneutral plasma of Ca^ in a linear 
quadrupole trap, the limit radius is given by the following relation where ujx/^ir is in 
MHz: 

R = \ —^r- X 1-31 X 10"^ (40) 
V 2L Ux/2Tr 

Another issue concerns the phase transition to crystals. In finite size samples like 
we are dealing with, the transition from the liquid to the crystal state is not abrupt. 
Molecular dynamics simulations [12] and experiments [33] have shown that it takes 
place for values of the coupling parameter F between 150 and 200. As the ions arrange 
themselves into concentric spheroidal shells, the diffusion inside a shell persists even for 
r as big as 300 [I2] . To study experimentally such phase transitions in a finite system, 
it is then important to make sure that F as high as a few hundred can be reached. For a 
laser cooled dense sample, F reaches a maximum limit Fc depending on the limit density 
He { a ^ (A/^TrncY^^) and the temperature reached by the cooling: 

This relation shows that the maximum F that can be reached in an experiment depends 
on the temperature of the sample, and the harmonic pseudopotential rnu^. As an 
example, for calcium ions and for ujx/2'n' in MHz, 

F, = 4.8 ^ ' . (42) 



Several experiments [351 EH] with laser cooled ions in linear radiofrequency traps have 
shown that because of heating processes like radiofrequency heating, the Doppler limit 
is not reached with large samples and that an optimistic estimation for the limit 
temperature is 10 mK. In this condition, one has to make sure to produce a high enough 
ujx to be able to explore a wide range of F. 

All the behaviours we have just described for quadrupole linear trap can not 
be extrapolated to multipole traps because the cold fluid limit does not result in a 
uniform density but rather in a profile density increasing with the radius (see Eq 1271). 
Nevertheless, we can reasonably assume a thermal equilibrium and a Gibbs distribution 
for multipole also. Indeed, the experimental results of the evaporation of ions from a 
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linear 22-pole [7] are consistent with a Boltzmann statistics to characterise the velocity 
distribution. 

3.2. In a linear multipole 

In this section, we use the method developed for the quadrupole linear trap in the 
more general case of a multipole linear trap. The differential equation governing the 
logarithmic density can be also written for a multipole potential, the major difference 
is that now the Laplacian of the pseudopotential depends on r: 

(t^trap{r, z) = ^{-] - -uy + -yy. (43) 



32^fc \rj Aq ' 2q 

With the reduced parameters (p, defined above, the general differential equation 
(!26|) for a multipole is 



no 32£krQ 

= exp[vl/(p,0]- V'^-^ (45) 

In practice, the numerical integration depends only on a = {eoVQ{2k — 
2)^A^~'')/(rao32£fero'^~^). The use of reduced parameter must not hide the limit of 
this method which refer the density profile to the central density uq- This is certainly a 
limitation here as, from molecular dynamics simulations [31], an empty core is expected 
in the cold limit. We see further in the text how to interpret results obtained from this 
model. As the integration of does not give access directly to the central density uq, 
to scale the size of the sample one has to use the same trick as in section 13. ![ which is 
to find a combination of a and Ad independent of uq. From the definition of A^, it is 
easy to see that a scales like nl~^, so the good combination is 



l2 /QOL rrc „2fc-2\ l/e^-l) 



A^ 32kBTSkr'o 



\{2k-2)^q^Vi) ' ^^^^ 

The idea of the numerical integration is the same as in section 13.11 First we assume 
a translation symmetry along Oz and do the integration along p only. Second, as the 
total number of ions per unit length is always defined by equation (IMl) . the three coupled 
parameters are now T, a and N/2L. The smaller a is, the higher is the number of ions 
for a given temperature. Then, the real size of the cloud is defined by XoPmax- Figure [3] 
shows the density profile of the same cloud in the same octupole trap but for different 
temperature. For a matter of numerical integration, only profiles for clouds as big as 
1.6 X 10^ ions/mm at 5 K could be calculated. For denser or colder sample, the profile is 
so steep that the adaptative step procedure I use for integration reached the lowest step 
allowed by the subroutine [37]. This limitation explaines why high reduced parameter 
p can not be reached by the simulation, contrary to the quadrupole case. The main 
difference from the quadrupole geometry is visible for the case of lower temperature 
(e,f) where the density reaches a maximum more than 20 times higher than the central 
density. This illustrates also the limitation we mentionned above for the cold fluid model. 
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in comparison with molecular dynamics simulations [21] • Choosing the density in the 
center of the trap uq as a reference for the density profile prevents one to find an empty 
core but the high ratio between the maximum density and no can be interpreted as a 
demonstration of this phenomena. This effect increases with the order of the multipole 




Figure 3. Density profiles versus tlie reduced radius p (a,c,e) or tlie radius r (b,d,f) 
for a prolate cloud with 1.6 x 10'' ions per mm, at different temperatures (a,b: T=10 
000 K, a = 19000, c,d: T=300 K, a = 1.8, e,f: T=5 K, a = 0.13). The scaling factor 
A D between p and r is calculated for Ca+ ions in an octopole with ro = 1 cm, Vo = 800 
V and ri/27r = 10 MHz. Notice that the x-axis of figure f has a different scale from b 
and d. 

as can be seen on figure H] which shows a comparison between the density profile in a 
8-pole and a 12-pole, for the same number of ions per length and the same trapping 
parameters. In the higher order geometry, the radial size of the cloud is 20% bigger 
than in the lower order geometry and the ratio of the maximum density over the central 
density is two times bigger. For the validity of our cold fluid model, it is worth noting 
that, by analogy with the quadrupole case, we defined the Debye length A^i relatively 
to the central density, which results into high values for X^, compared to the radial size 
of the cloud (see figure [3] for values). Nevertheless, the relevance of our model depends 
on the local Debye length Xd{p) compared to the size of the sample. We can see on 
figure m that, close to the edge of the cloud, this length Xd{p) is a lot smaller than Xd, 
as the local density is far higher than uq. Consequently, using the cold fluid model to 
scale such a cloud seems valid in the cold and dense limit. 

As for a quadrupole, the cold fluid limit allows one to estimate the minimum size of 
a cloud by assuming the T ^ limit. To prevent divergence of the density at this limit, 
exp[\l/(p)] — ap^^~^ —>■ 0, with our previous notation. The maximum reduced radius 
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1 2 3 4 5 6 

Km) xlO"^ 



Figure 4. Density profile of a calcium cloud of 1.6 x lO** ions/mm at 5 K in two traps 
of different order but with same parameters as for figure[3]: ro = 1 cm, Vq = 800 V and 
f]/27r = 10 MHz. Dashed line: fc = 4, a = 0.13, Xd = 0.86 mm, uq = 3.2 x 10^" m^^. 
Solid line: k = 6,a = 0.0056, Xd = 1-4 mm, m = 1.2 x 10^'' m^^. 



visited by the ions is then determined by the hnear density N/2L: 

The dependance of the cold fluid limit radius Rm on the trapping parameters can then 
be expressed as 

Rm = ro(^ . (48) 

To put values into this equation, let's choose the case of calcium ions in an octopole 
{k = 4) and express Q/2tt in MHz, then 

M) [ Vo 

As an example, for Vq = 400 V, Q/2'k = 1 MHz and r^ = 1 cm, the radial size of a 
cylindrical sample of N/2L = 4.2 x 10^ ions/mm is Rm = 2.4 mm in the cold fluid limit. 
As a comparison, at 300K, the numerical integration gives R = 3.8 mm. 

Also like in quadrupole traps, the set of trapping parameters is constrained by 
stability criterion for the trajectories of the ions. But contrary to quadrupole traps, 
there is no absolute criterion but one can use the local adiabaticity criterion of Teloy 
and Gerlich introduced in 12.11 For a perfect multipole like we are assuming, 

Vad = -^TT^ - . (50) 



Rm={7rr) I \/ " 1 xo.45. (49) 



2k£k \ro^ 

To be sure to remain in the adiabatic regime, we want to keep r]ad < Vitm and then the 
radial size R < r^^^ with 

/ nj,C \l/(fc-2) 
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To have an idea of how this criterion is respected, it can be useful to compare, for a 
given set of trapping parameters, the smallest achievable radial size Rm and the highest 
allowed radius for adiabatic trajectories 



For Calcium ions in an octopole, this becomes 

— 1— = 5 — TTTT X 12.9. (53) 

As an example, for the same parameters as above (Vq = 400 V, i7/27r = 1 MHz, tq = 1 
cm) and for the same number of ions per unit length N/2L = 4.2 x 10^ ions/mm, 
Rm/'^max — ^-"^^ ^^"^ adiabaticity criterion is obeyed over all the cold sample (we 
used r]iim = 0.3 here). 

By looking at equations (I46f48ll52l) . one can notice that if the ratio fi/Vo is 
conserved, the absolute radial size is also conserved but not the ratio Rm/r'^ax- 
practice, by substituting Vl/Vq by 2n/2Vo, Rm/i"max is reduced by a factor \/2. One 
can use this difference in behaviour to scale ones trap and make sure that the desired 
sample fits inside the adiabatic volume. 



4. Conclusion 



In this tutorial we have described the dynamics of charged particles in a radiofrequency 
trap in a very general manner to point out the differences between the dynamics in a 
quadrupole and in a multipole trap. When dense samples are trapped, the dynamics is 
modified by the Coulomb repulsion between ions. To take into account this repulsion 
into the equilibrium state of the cloud, we use a method that models the ion cloud as a 
cold fluid. In the case of prolate clouds, it allows one to scale the size of the samples with 
the trapping parameters and the number of trapped ions, for different linear geometries 
of trap. We think this can be useful to build an experiment where a large number of 
ions needs to be trapped. 
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